This is the third paper of a series devoted to study the stel- 
lar content of early-type galaxies. The goal of the series is 
to set constraints on the evolutionary status of these objects. 
We present line-strength gradients for 22 spectral indices mea- 
sured in a sample of 82 early-type galaxies in different environ- 
ments, including the high-density core of the Coma cluster, the 
Virgo cluster, poor groups, and isolated field galaxies. Using 
new evolutionary population synthesis models we derive age 
and metallicity gradients, and compare the mean values with 
the predictions of different galaxy formation models. We ex- 
plore the behaviour of individual chemical species by deriv- 
ing the metallicity gradient with different indicators. We find 
that the strength of the metallicity gradient inferred from stel- 
^O ■ lar population models depends on the specific Lick index em- 
^D ployed. In particular, metallicity gradients obtained with CN2 
^— ^ and Ca4668 combined with HjS are steeper than when measured 
' using Ca4227 or Fe4383. The correlation of the metallicity gra- 
^ dients with other parameters also depends on the specific index 
.^ employed. If the metallicity gradient is obtained using CN2 and 
Mgb then it correlates with the central age of the galaxies. On 
t^ the contrary, if Fe4383 or Ca4227 are used, the metallicity gra- 
^^ dient correlates with the velocity dispersion gradient. This may 
,__! ' suggests that several mechanism have helped to set the age and 
J> metallicity gradients in early-type galaxies. While we do not 
T-H _ find any correlation between the metallicity gradient and the 
C ■ central velocity dispersion for galaxies in low-density environ- 
2_? ' ments, we find a marginal correlation between the metallicity 
^-v gradient and the mass for galaxies in the centre of the Coma 
\^ ' cluster. We also find a trend for which galaxies in denser en- 
^^ vironments show a steeper metallicity gradient than galaxies 
(— I , in less dense environments. We interpret these results in light 
Q_i' of the different mechanisms proposed to explain the observed 
' ' changes between galaxies as a function of environment. 
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Abstract. 



1. Introduction 

This is the third paper in a series devoted to studying the prop- 
erties of the stellar populations in early-type galaxies as a func- 
tion of their local environment. In Sanchez-Blazquez et al. 
(2006a, hereafter Paper I) we presented central Lick/IDS in- 
dex measurements for a sample of 98 early-type galaxies and 
described their relation to the velocity dispersion. In Sanchez- 
Blazquez et al. (2006b, hereafter Paper II) we compared the in- 
dices with the stellar population synthesis models of Vazdekis 
et al. (2006, in preparation; hereafter V06) to derive ages and 
metallicities. These models are an improved version of those 
described by Vazdekis et al. (2003) strengthened by the ad- 
dition of a new stellar library (Sanchez-Blazquez et al. 2006, 
in preparation). The synthetic spectra cover a spectral range of 
/1/13500-7500 A at resolution 2.3 A. The previous papers con- 
centrated on the analysis of the central regions of the galaxies. 
In this paper we investigate the behaviour of these properties as 
a function of galactocentric radius. 

It is well known that early-type galaxies show a variation 
of their stellar population properties with radius. The first ev- 
idence of this phenomenon was observed in the colours of 
bulges and elliptical galaxies (see the review by Kormendy & 
Djorgovski 1989) and indicated that the central regions tend 
to be redder than the outer parts. Other authors have measured 
optical surface brightness profiles for large sets of local ellip- 
ticals (e.g. Franx, Illingworth & Heckman 1989; Peletier et al. 
1990), reaching the same conclusion. The first study of gradi- 
ents in the strength of absorption features was performed by 
McClure (1969). In a sample of 7 galaxies, McClure found that 
the C(41 -42) index (a measure of the CN/14216 band strength) 
was stronger in the centres of the galaxies than at a distance 
of 1-1.5 kpc, which he interpreted as a metallicity difference 



between the two regions. Subsequent work in the field has 
explored the strength of absorption feature gradients using a 
broad range of line indices (Spinrad et al. 1971; Spinrad, Smith 
& Taylor 1972; Welch & Forrester 1972; Joly & Andrillat 
1973; Oke & Schwarzschild 1975; Cohen 1979; Efstathiou & 
Gorgas 1985; Couture & Hardy 1988; Peletier 1989; Thomsen 
& Baum 1989; Gorgas, Efstathiou & Aragon-Salamanca 1990; 
Boroson & Thompson 1991; Bender & Surma 1992; Davidge 
1992; Davies, Sadler & Peletier 1993; Carollo, Danziger & 
Buson 1993; Gonzalez 1993; Fisher, Franx & Illingworth 1995; 
Gorgas et al. 1997; Cardiel, Gorgas & Aragon-Salamanca 
1998a; Mehlert et al. 2003). In general, these studies indicate 
the existence of intense gradients in CN/13883 and CN/14216, 
less pronounced gradients in the Mg/I5176, G band, NaD, Ca 
H&K features, and in some lines of Fe i, and a weak or null 
gradient in H/?, MgH, TiO, Ca i and in the calcium triplet in 
the near-infrared. Most of these studies have suggested that the 
existence of gradients in the metallic spectral features is a con- 
sequence of a decreasing metallicity with increasing galacto- 
centric radius (e.g. Mc Clure 1969; Cohen 1979; Davies et al. 
1993; Kobayashi & Arimoto 1999; Mehlert et al. 2003). Some, 
however, have argued that, apart from a variation of metallicity 
with radius, there is also a radial variation in the luminosity- 
weighted mean age of the stellar populations, with the cen- 
tral regions being younger than the outer regions (Gorgas et al. 
1990; Munn 1992; Gonzalez 1993; Gonzalez & Gorgas 1996). 

How the physical properties of galaxies vary with radius 
can prove invaluable for constraining the processes of galaxy 
formation and evolution. For example, metallicity gradients are 
a measure of the quantity, velocity, and duration of gas dissipa- 
tion. Likewise, age and metallicity gradients contain informa- 
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tion concerning the relative importance of interactions during 
galaxy formation. 

In the classical models of monolithic collapse (Eggen et al. 
1962; Larson 1974a; Carlberg 1984; Arimoto & Yoshii 1987; 
Gibson 1997), stars form in essentially all regions during the 
collapse and remain in their orbits with little inward migra- 
tion, whereas the gas dissipates inwards, being continuously 
enriched by the evolving stars. In this way, the stars formed 
at the centres of galaxies are predicted to be more metal-rich 
than those born in the outer regions. Supernova-driven galac- 
tic winds (Mathews & Baker 1971; Larson 1974b; Arimoto & 
Yoshii 1987; Gibson 1997), initiated when the energy injected 
into the interstellar medium (ISM) by supernovae matches that 
of its binding energy, act to evacuate the galaxy of gas, thereby 
eliminating the fuel necessary for star formation. The exter- 
nal parts of the galaxy (with a shallower potential well) de- 
velop winds before the central regions, where the star for- 
mation and, therefore, the chemical enrichment continue for 
longer. The monolithic collapse models, therefore, predict very 
steep metallicity gradients as both processes - the dissipation 
of gas toward the central parts of the galaxy and the different 
timescales for the occurrence of the galactic winds - act in the 
direction of steepening any nascent metallicity gradient. ' 

Simulations of galaxy mergers within the concordant hi- 
erarchical clustering cold dark matter framework (e.g. Cole 
et al. 1994; Baugh, Cole & Frenk 1996; Kauffmann 1996; 
Kauffmann & Chariot 1998) offer somewhat contradictory pre- 
dictions as to the radial variation of stellar properties in early- 
type galaxies - while some (e.g. White 1980; Bekki & Shioya 
1999) suggest mergers lead of a flattening of metallicity gradi- 
ents, others (e.g. van Albada 1982) argue that the gradients are 
affected only moderately by the mergers, as the violent relax- 
ation preserves the position of the stars in the local potential. 
This apparent dichotomy between the simulations is driven in 
part by the sensitivity of the outcome to the fraction of gaseous 
versus stellar mass present in the progenitor galaxies. Broadly 
speaking, the predicted gradients are steeper if the progenitor 
galaxies have a large fraction of their pre-merger baryonic mass 
in the form of gas. Furthermore, numerical simulations suggest 
that during the merger a significant fraction of this gas migrates 
inward toward the central regions of the merging galaxies, re- 
sulting in increased central star formation (Barnes & Hernquist 
1991). Mihos & Hernquist (1994) showed that the observed 
gradients in elliptical galaxies may be a consequence of the 
occurrence of secondary bursts of star formation triggered by 
these mergers. 

One of the keys to dermining the physical underlying the 
formation and evolution of galaxies is to study of the relations 
between the gradients and other fundamental (global) prop- 
erties of galaxies. For instance, dissipational collapse models 
predict a strong positive correlation between metallicity gradi- 
ent and galactic mass. The empirical evidence for such puta- 



' Galactic winds can produce metallicity gradients without invok- 
ing dissipation (e.g. Franx & Illingworth 1990; Martinelli, Matteucci 
& Colafrancesco 1998), with the local metallicity coupled directly to 
the local potential well depth and essentially independent of the galaxy 
formation collapse physics. 



tive correlations remains contentious - e.g. Gorgas et al. (1990, 
G90 hereafter) did not find any relation between Mg2 gradi- 
ent and the rotation or total luminosity of a sample of early- 
type galaxies, although they found some evidence of a positive 
correlation between Mg2 gradient and central velocity disper- 
sion (a probability of 95% in a non-parametric Spearman test). 
Conversely, Franx & Illingworth (1990) found a correlation 
between colour gradient and local escape velocity (later con- 
firmed by other authors, as Davies et al. 1993), arguing that the 
aforementioned galactic winds were the dominant mechanism 
controlling the metal content of early-type galaxies. Davidge 
(1992) analysed 12 bright ellipticals, searching for correlations 
between Mg2 gradients and the central velocity dispersion (o-q), 
the total luminosity, the shape of the isophotes, the fine struc- 
ture parameter (Schweizer et al. 1990), and the anisotropy pa- 
rameter (v/cr)*, which defines the degree of rotation in a galaxy 
(Binney 1978). Davidge found weak correlations between the 
Mg2 gradient and (v/cr)* and ctq, and an absence of correlation 
with any of the other parameters. CaroUo et al. (1993) carried 
out an exhaustive study of the gradients of several spectral fea- 
tures (Mgi, Mg2, NaD, TiOi, Ti02, and Fe5270) in a sample 
of 42 galaxies. Carollo et al. found a tendency for the slope 
of the gradients in the Mg2 index to increase with the mass of 
the galaxy, but only for galaxies with masses below lO" Mq. 
Carollo & Danziger (1994) studied the line-strength gradients 
in Mg2 and (Fe) for five early-type galaxies, confirming the 
dependency between metallicity and local potential well depth 
as a function of galactocentric radius, albeit with significant 
scatter Gonzalez & Gorgas (1996) surveyed the Mg2 gradi- 
ent literature available at the time, discovering that galaxies 
with steeper Mg2 gradients also possess stronger central Mg2. 
They proposed a scenario in which star formation episodes in 
the centre of the galaxies are responsible for the correlation, 
pointing out that its existence implies that the global mass- 
metallicity relation is flatter than the mass-metallicity relation 
inferred from the central values. Recently, Mehlert et al. (2003) 
found a correlation between the gradients of some spectral fea- 
tures and the velocity dispersion gradient (which can be consid- 
ered a measure of the potential well depth gradient). In contrast 
with Gonzalez & Gorgas, Mehlert et al. did not find a correla- 
tion between the line-strength index gradients and the central 
values, or between the index gradients and the central velocity 
dispersion. 

The results to date are, to some extent, contradictory and 
therefore incapable of providing unequivocal support to any 
particular galaxy formation scenario. This contradictory nature 
is driven, in part, by the requisite high signal-to-noise data and 
associated care in data reduction necessary to extract reliable 
gradients. Further, very few absorption features (mainly Mg2) 
have been systematically explored with a suitably large sam- 
ple of ellipticals. Our work has been designed specifically to 
address these shortcomings which have plagued the interpreta- 
tion of the extant data, using a sample of 82 galaxies populating 
a range of local environmental conditions. 

In Section 2, we describe the measurement of the line- 
strength gradients within our dataset. In Section 3, we com- 
pare these gradients with synthesis models to derive explicit 
age and metallicity gradients. Sections 4 to 7 analyse the stel- 
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lar population gradients for a sample drawn from low den- 
sity environments. In particular, Section 4 discusses the mean 
age and metallicity gradients, while in Section 5 we explore 
the existence of possible variations of the chemical abundance 
ratios with radius. In Section 6, we explore the existence of 
putative correlations between metallicity gradient and other 
galactic parameters. In Section 7, we use the stellar popula- 
tion global parameters, derived with the help of the gradients, 
to study whether the trends found in Papers I and II hold only 
for the central parts, or can be extended to the entire galaxy. 
Section 8 studies the mean gradients in the sample of galaxies 
drawn from the high-density core of the Coma cluster. Finally, 
in Section 9, we summarise our findings and conclusions. 

2. Measurement of gradients 

The parent sample from which galaxies are drawn consists of 
98 early-type systems, of which 37 belong to the Coma cluster 
(high-density environment galaxies, hereafter HDEGs), while 
the rest are drawn from the field, small groups, and the Virgo 
cluster (low-density environment galaxies, hereafter LDEGs). 
Details of the sample, the observations, and associated data 
reduction, can be found in Paper I. From each fully reduced 
galaxy frame, a final frame was created by extracting spectra 
along the slit, binning in the spatial direction to guarantee a 
minimum signal-to-noise ratio per A (S/N) of 20 in the spec- 
tral region of the H/3 index, ensuring a maximum relative error 
for this index (assuming a typical value of 1 .5 A) of 20% (see 
Cardiel et al. 1998b). Those galaxies which did not have a min- 
imum of 4 spectra along the radius meeting this criterion were 
eliminated from the final sample, which left a total of 82 galax- 
ies of sufficient quality to measure gradients - 21 HDEGs and 
61 LDEGs. 

An equivalent radius ?(i) was assigned to each binned spec- 
trum as 

j.fj.i{r)rdr 



m 



(1) 



j.Mi{r)dr ' 

where 

Mr) = k lO*'--. (2) 

The value of /?, for each spectrum was calculated by fitting the 
number of counts in each row of the detector (including all the 
rows in each spectrum plus the neighbours) to the relation 

logN ^bir + c, (3) 

where A^ is the number of counts, r the radius in arc- 
seconds, and c a constant. We projected the gradients over 
the major axis, multiplying the mean radius by a factor 

/ - Jcos^ ^ + (f ^™'^ ^)' where 6 is the difference between the 
major axis and the position angle of the slit, and a and b rep- 
resent the major and minor axis of the galaxy respectively. The 
a, b and 9 values for each galaxy can be found in Table 2 of 
Paper I. 

Radial velocities and velocity dispersions were measured 
for each spectrum as a function of galactocentric radius as de- 
scribed in Paper I. The rotation curve and the velocity disper- 
sion profile were used to measure the Lick/IDS indices at the 



resolution defined by this system, and to correct them for veloc- 
ity broadening (again, as described in Paper I). Instead of using 
the individual values of the velocity dispersion for each spec- 
trum measured at each radial "bin", we fitted a smooth curve 
to the cr profile. This minimises the errors in the derived ve- 
locity dispersion, due to the lower signal-to-noise ratios in the 
external parts of the galaxies. Residual [OIII] emission pro- 
files were also estimated as described in Paper I. These val- 
ues were used to correct the iUfi index for emission using the 
the relation derived by Trager et al. (2000a) - AH8=0.6 [OIII], 
where AHjS is the correction to the index, and [OIII] the equiv- 
alent width of the [OIII]/15007 emission line. The Mgb and 
Fe5015 indices were also corrected from emission, as in the 
central spectra. The final kinematic, line-strength and emission 
profiles, together with their associated errors, are presented in 
Sanchez-Blazquez (2004). 

We measured 19 Lick/IDS indices defined by Trager et al. 
(1998) for all the galaxies except for those observed in Run 3, 
for which only 15 indices could be measured. With the aim of 
reducing the random errors in the external parts of the galaxies, 
we have defined the following composite indices: 



HA = 



Balmer 



HyA + mA 

2 
HjA + H^A + HyS 



(4) 
(5) 



To obtain the final values for the line strength indices we 
applied the same corrections as for the central indices (see 
Paper I). Finally, all indices were transformed into magnitudes, 
as described in Paper I. 

Although the behaviour of the indices as a function of 
galactocentric radius is not perfectly linear for all 82 galax- 
ies, with the aim of obtaining a quantitative measurement of 
their values, we performed a linear fit, weighted by the errors 
over all the indices, to the relation I' - c + d\og — , where 
/' represents the index expressed in magnitudes, a^ is the ef- 
fective radius projected over the major axis, i.e. Ue - r^ffj J^ 
(reff being the eff'ective radius), and d the index gradient, here- 
after denoted by grad. This projection can be done under the 
assumption that the contours of constant /' coincide with the 
isophotes, which has been confirmed by several earlier studies 
(e.g. Davies et al. 1993). To avoid the uncertainties and radial 
flattening due to seeing eff'ects, the data within the central area 
corresponding to the seeing were excluded from the fit. The 
seeing was estimated directly from the observations on the 4 
diff'erentruns. Fig.[nshows some examples of the line-strength 
gradients for the galaxy NGC 2832. 

In order to check the assumption of linearity for the gra- 
dients we have compared the central values derived in Paper I 
(integrated within an equivalent aperture of 4" at redshift z - 
0.016 - i.e., ~1.3 kpc), with the expected values for this aper- 
ture inferred from the gradients. Fig.|2lshows this comparison 
for all the indices. Within each panel, the mean offset (A) and 
the root mean square of the deviations (cr) are indicated. In gen- 
eral, the agreement is very good - the offsets are not statistically 
significant for any of the indices and the scatter is compatible 
with the errors, lending support to the adopted approximation 
of a linear fit. 
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Fig. 1. Variation of the indices CN2, C4668 and H/3 as a function of galactocentric radius for the galaxy NGC 2832. The upper 
axes show the radial distance from the centre of the galaxy in arcseconds. The solid line represents a linear fit to the data, as 
described in the text. The data within 0".9 of the galaxy's centre were excluded from the fit to minimise the effects of seeing on 
the fitted gradients. 
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Fig. 2. Comparison of the Lick indices measured in the central spectra (integrated within an equivalent aperture of 4" at redshift 
Z = 0.016 - i.e., --1.3 kpc) and the indices obtained from the gradients integrating over the same aperture; cr indicates the root 
mean square and A the mean offset between both measurements. The solid line indicates the 1 : 1 correspondence. In all cases, the 
mean offsets are not statistically significant. 
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Table n] available in the electronic edition and at the web 
site^, lists the gradients in all the line-strength indices measured 
for each of the 82 galaxies in our sample. A portion of this table 
is reproduced here to demonstrate its content and structure. 



3. Gradients in age and metallicity 

Figure|3lshows the line-strength gradients in Fe4383, Mgb and 
HyS, presented in index-index diagrams that combine both the 
central values and the line-strengths at one effective radius. 
Each central measurement (solid circle) is connected to its cor- 
responding measurement in the outer region. Overplotted are 
the stellar population models of V06. As can be seen, most of 
the lines tend to be nearly parallel (although, admittedly, not 
entirely) to the iso-age lines taken from the theoretical grid, 
which indicates that the gradients are mostly due to variations 
in metallicity. However, there are several galaxies which show 
significant variations in both age and metallicity with radius. 
We return to this point towards the end of this Section. 

To quantify this behaviour, we have transformed the line- 
strength gradients into age and metallicity gradients using 
the model predictions of V06, who make use of the code of 
Vazdekis (1999) updated with a new and improved stellar li- 
brary (MILES, Sanchez-Blazquez et al. 2006; these models 
will be presented in the forthcoming paper V06). 

With an aim to minimising the effects of low signal-to-noise 
data, we have designed a method that makes simultaneous use 
of 10 different indices. The selected indices are: H6a, CN2, 
Ca4227, G4300, Hja, Fe4383, H/3, Fe5015, HA, and Balmer. 
We checked that the results did not depend on the particular 
choice of indices. None of the Mg indices were employed, as 
they could not been measured in all the galaxies from Run 3 
due to the wavelength converage of the spectra. In order to de- 
rive single estimates of age and metallicity gradients, we have 
followed a procedure which is divided into two steps: (i) ob- 
taining individual errors for the age and metallicity gradients 
for all the galaxies in each index-index diagram, and (ii) deriv- 
ing the final values of the age and metallicity gradients, along 
with their associated errors. 




Mgb 




Fe4383 

Fig. 3. Index-index diagrams using the models of V06. Solid 
lines represent the predictions for populations of constant age 
(top to the bottom, 1.41, 2.00, 2.82, 3.98, 5.62, 7.98, 11.22, 
15.85 Gyr), while dashed lines show populations of constant 
metalhcity {left to right, [M/H]=-0.68, -0.38, 0.0, -hO.2). The 
solid circles in the figure indicate the position of the centre of 
the galaxies, while the lines connect with the measurements at 
one effective radius. 



3.1. Obtaining ttie errors in tine age and metallicity 
gradients 

We first calculated the slopes of the constant age and metallic- 
ity lines in the models. The slopes of these Unes change with 
the absolute values of age and metallicity but, in the region in 
which the galaxies are located, the values are almost constant. 
The derived slopes define a new diagram that we call "Aindex- 
Aindex" - in the associated Fig. 0] the dashed line represents 
the expected gradients in the indices if these were due, exclu- 
sively, to a variation in age (assuming a constant solar metal- 
licity). On the other hand, the solid lines show the predicted 
gradients if the only parameter changing with radius was the 
metallicity (assuming a constant age of 10 Gyr). Fig. |3 also 
shows the gradients in H/3 and Fe4383 measured in the galaxy 
NGC 4842A. To derive the age and metallicity gradients, we 



' ^http : //www ■ ucm . es/info/Astrof /users/pat/pat . html | 



project the position of the galaxy over the lines Aage=0 and 
A[M/H]=0, as shown by the dotted lines in the figure, and in- 
terpolate the projection over these lines. 

For each galaxy, and in every Aindex-Aindex diagram, we 
performed lO'* Monte Carlo simulations in which each point 
was perturbed with our errors, following a Gaussian probabil- 
ity distribution. For each simulation, an age and metallicity gra- 
dient was obtained. This process was repeated with all possi- 
ble paired combinations of indices. As an illustration. Fig. |5] 
shows the values of age and metallicity gradients derived for 
NGC 4842A, using the AH^-AMgb and AMgb-AFe4383 di- 
agrams. As can be seen, due to the non-orthogonality of the 
index-index diagrams, there is an artificial anti-correlation be- 
tween the age and the metallicity gradients. To obtain the error 
in the age and metallicity gradient for each galaxy, in every di- 
agram, we projected these ellipses over the x- and y-axes and 
calculated the standard deviation of the resultant Gaussian. We 
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Table 1. Velocity dispersion and line-strength gradients for the sample of 82 galaxies. A portion of the table is shown for guidance 
regarding its form and content. For each galaxy and index, the first line shows the measured gradient, while the second row lists 
its corresponding formal error The full table is available in the electronic edition of this paper. 
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Fig. 4. Variations in the H/3 and Fe4383 line-strength indices 
expected through changes in metallicity at a constant age of 
10 Gyr (solid line), and by variations in the age at constant solar 
metallicity (dashed line). The point represents the values of the 
gradients for the galaxy NGC 4842A with their associated error 
bars. The dotted lines indicate the projections of the values over 
the lines of constant age and constant metallicity. 




denote the typical errors by cr(grad[age]), and cr(grad[M/H]),, 
where / refers to each index-index diagram. 



3.2. Obtaining the age and metallicity gradients 



Once we have obtained the typical errors associated with the in- 
ferred age and metallicity gradients for each index-index com- 
bination, we again performed (for each galaxy) 10"* Monte 
Carlo simulations, into which Gaussian noise was added to 
each line-strength index. Simulations were performed for each 
index independently, as opposed to each iindex-index diagram 
(note that once an index value is simulated, its particular value 
is fed simultaneously to all the possible diagrams taking into 
account the fact that not all of the different diagrams are inde- 
pendent). 

We then measured the age and metallicity gradients of all 
the values obtained in all the simulations j, using all the dif- grad[age], - - 
ferent diagrams /. The mean age and metallicity gradient of 
simulation j was obtained as an average of the age and metal- 



5 

grad(age) (Fe4383-Mgb) 

Fig. 5. Age and metallicity gradients obtained for the 10"* 
Monte Carlo simulations performed in the diagrams AH/?- 
AMgb (top panel) and AMgb-AFe4383 (bottom panel) for 
NGC 4842A. The correlation of the errors is larger when the 
age and metallicity are inferred from the latter. 



licity gradients obtained for this particular simulation using all 
the different diagrams i, namely 



grad[age],.^. 
^ cr(grad[age])2 



"com 



"com 

Z 

1=1 



1 



(6) 



cr(grad[age]); 
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and 



grad[M/H],- = f 



grad[M/H],;,. 
.^j .(grad[M/H])2 



y grad[ 
^ cr(grad 



(7) 



Z 



i cr(grad[M/H])f 



where grad[age],j- represents the age gradient obtained with the 
f^ diagram in the simulation j, cr(grad[age]), the uncertainty of 
the age gradient in /-th diagram (obtained in the first step of 
the process), grad[M/H],-j the metallicity gradient in simula- 
tion /■ (using the «* diagram), and cr(grad[M/H]), the error in 
the metallicity gradient in this diagram; ncom indicates the to- 
tal number of diagrams, which is the number of possible index 
pairs (ncom=45). 

Finally, the age and metallicity gradients for each galaxy 
were calculated as the mean values of all the simulations, i.e. 



grad[age]fi„^i = ^ 



grad[age^] 



and 



grad[M/H]fi„,, = \ i, 

■^ «siin 



(8) 



(9) 



being «sin, the number of simulations performed for each line- 
strength index gradient. The final errors in the age and metal- 
licity gradients were obtained as the standard deviations of the 
values obtained for each simulation, that is 



"siin 

^ (grad[age]^. - grad[age]g„^,) 



0"(age) - \ 

and 



;=i 



'^sim ^ 



(10) 



0-([M/H]) - \^ 



2 (grad[M/H]^. - grad[M/H]fi„„)' 



j=i 



(11) 



This process was performed for all 82 galaxies in our sam- 
ple. The final values for the age and metallicity gradients ob- 
tained in this way, together with the associated errors, are listed 
in Table 121 

Figure |6l shows three difi'erent Aindex-Aindex diagrams in 
which we have over-plotted the galaxies in low-density envi- 
ronments. The lines of constant age and metallicity divide these 
diagrams into four distinct regimes. These regimes indicate the 
position of the galaxies with positive/negative differences in 
age and metallicity between the central and the external parts, 
as labeled in the panels. As can be seen, most of the galax- 
ies are situated in the regime consistent with their centres be- 
ing younger and more metal-rich than their outer regions. The 
bulk of the galaxies, however, are further away from the line 
A[M/H]=0 (dashed line) than from the line Aage=0 (solid line), 
suggesting that the gradients in the line-strength indices are due 
(mostly) to variations of metallicity with radius. 



center less metallic end older 



center, more ^ 

metollic and older ^ 




ceriter, I 
metollic 
end younger. 



center more metollic and younger 



-0.2 



-0.1 



0.1 



AMgb 



X 




center more metollic and younger 



_1_ 



-0.2 



-0.1 

ACN2 



0.1 



center less metallic end older 



ceriter more , ,^ , 
metallic and older 




center more metallic and younger 
_i I I I . ^ (. . . 



-0.1 



0.1 



AFe4383 



Fig. 6. Aindex-Aindex diagrams combining the Hj6 gradient 
with Mgb, CN2, and Fe4383. Solid lines indicate the slope of 
the predictions at constant age (10 Gyr); the dashed line shows 
the slope of the simple stellar population at solar metallicity. 
The solid circles represent the values of the line-strength gra- 
dients for the sample of galaxies in low-density environments. 



4. Mean gradients 

The predicted mean values of the age and metallicity gradi- 
ents are dependent upon the merging history of the galaxies. 
In general, dissipative processes tend to steepen the metallic- 
ity gradient while major mergers are expected to dilute it. In 
the simulations of Kobayashi (2004), the typical gradients for 
non-merger and merger galaxies are A[Fe/H]/Alog r ~ -0.45, 

-0.38, A[0/H]Alogr 0.25, -0.24, and A log Z/ log r ~ 

-0.30, -0.24, respectively. We now analyse the mean age and 
metalhcity gradients for the sample of LDEGs. Table |3l list 
these values, together with relevant associated statistics (see 
below). 
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Table 2. Age and metallicity gradients derived with the method of Sec. |3 (second and third columns), and metallicity gradients 
derived in several index-index diagrams using different metallicity indicators (as indicated in the column headers) combined with 
HyS. The second row for each galaxy indicates the associated errors in the gradients. 









CN2 


C4668 


Fe4383 


Mgb 


Ca4227 


Galaxy 


grad(age) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 




±(T 


±ar 


±0- 


±cr 


±cr 


±(T 


±0- 


NGC 221 


0.0151 


-0.0273 


-0.170 


-0.246 


0.023 


-0.171 


0.082 




0.0180 


0.0226 


0.241 


0.139 


0.222 


0.259 


0.373 


NGC315 


0.0850 


-0.2745 


0.010 


-0.689 


-0.214 


-1.266 


-0.209 




0.0486 


0.0651 


0.112 


0.089 


0.095 


0.148 


0.151 


NGC 507 


-0.6291 


0.5756 


0.076 


0.451 


1.004 


0.453 


0.935 




0.1508 


0.2140 


0.309 


0.324 


0.354 


0.611 


0.664 


NGC 584 


0.0814 


-0.1937 


-0.673 


-0.510 


-0.341 


-0.414 


-0.294 




0.0868 


0.1185 


0.196 


0.114 


0.145 


0.242 


0.355 


NGC 636 


0.0779 


-0.3627 


-0.740 


-0.804 


-0.427 


-0.291 


0.117 




0.1148 


0.1592 


0.220 


0.200 


0.300 


0.274 


0.382 


NGC 821 


0.0563 


-0.8237 


-1.560 


-1.058 


-0.493 


-0.981 


-1.057 




0.1157 


0.9981 


0.351 


0.235 


0.293 


0.391 


0.478 


NGC 1453 


-0.2652 


-0.0141 


-0.170 


-0.246 


0.023 


-0.171 


0.082 




0.1036 


0.1259 


0.241 


0.139 


0.222 


0.259 


0.373 


NGC 1600 


0.2293 


-0.0732 


-0.765 


-0.521 


-0.592 


-0.561 


-0.897 




0.0504 


0.0708 


0.116 


0.072 


0.121 


0.216 


0.203 


NGC 1700 


0.0803 


-0.3432 


-0.750 


-0.693 


-0.096 


-0.384 


0.534 




0.1116 


0.1349 


0.186 


0.183 


0.237 


0.274 


0.508 


NGC 2300 


0.1531 


-0.2416 


-0.553 


-0.902 


-0.335 


-0.644 


-0.049 




0.1284 


0.1582 


0.270 


0.156 


0.252 


0.353 


0.508 


NGC 2329 


0.2413 


-0.3977 


-0.610 


-0.510 


-1.081 


1.243 


-0.014 




0.1731 


0.1840 


0.346 


0.284 


0.377 


0.884 


0.478 


NGC 2693 


-0.1917 


-0.0072 


-0.420 


-0.220 


0.049 


0.210 


-0.168 




0.0786 


0.0997 


0.116 


0.138 


0.206 


0.205 


0.265 


NGC 2694 


-0.1551 


-0.0214 


0.119 


-0.531 


0.270 


0.186 


-0.882 




0.2540 


0.2572 


1.146 


0.218 


0.467 


0.853 


1.057 


NGC 2778 


0.0644 


-0.3501 


-1.041 


-0.733 


-0.284 


-1.158 


-0.233 




0.0879 


0.1155 


0.186 


0.125 


0.178 


0.182 


0.249 


NGC 2832 


0.1436 


-0.3521 


-0.819 


-0.747 


-0.712 


4.801 


-1.200 




0.0696 


0.1046 


0.120 


0.102 


0.196 


1.205 


0.297 


NGC 3115 


0.2518 


-0.4864 


-1.121 


-0.691 


-0.416 


-0.692 


-0.367 




0.0501 


0.0845 


0.154 


0.070 


0.090 


0.103 


0.129 


NGC 3377 


0.0334 


-0.3600 


-1.118 


-0.960 


-0.498 


-1.143 


-0.229 




0.0679 


0.1147 


0.152 


0.103 


0.241 


0.214 


0.222 



4.1. Age gradients 

The mean age gradient for the sample of LDEGs is 
Alog(age)/logr=0.082 ±0.015. A f-test indicates that the 
probability of this value being different from zero by chance 
is < 0.1%. The positive gradient indicates that, on average, the 
centres of the galaxies are younger (or at least, they contain a 
percentage of younger stars which make the mean luminosity- 
weighted age lower) than the outer parts. The existence of a 
significant age gradient is difficult to explain within monolithic 
collapse scenarios, as the timescales for star formation are nec- 
essarily very short, but it does suggest that secondary episodes 
of star formation have occurred recently in the centres of the 
galaxies. 

Little work has been done to date on the derivation of age 
gradients within early-type galaxies. Indeed, most of the extant 
studies of gradients in the literature assume the absence of an 
age gradient in order to derive the metallicity gradient. Munn 



(1992) was the first author that studied, in a systematic way, 
the variation of age with the galactocentric radius. From a sam- 
ple of seven early-type galaxies, Munn combined the CN/13883 
and CN/14216 features with the D4000 index to compare with 
the prediction of stellar population models. Since the age cali- 
bration employed by Munn was necessarily restricted to the as- 
sumption of solar metallicity, his ability to quantify the inferred 
gradients was limited. He did however conclude that in order to 
explain the dispersion and trends in the CN-D4000 diagrams, 
a variation of at least two parameters with radius was necessary 
(as suggested already by G90), with age and metallicity being 
the obvious candidates. Using the stellar population models of 
Worthey (1994), Gonzalez (1993) found a mean variation in the 
age of his sample of galaxies of ~ 20% from the centre to the 
effective radius (he also found a variation of ~50% in metallic- 
ity). Fisher, Franx & Illingworth (1996) found in their sample 
of lenticulars that the centres were slightly younger than the ex- 
ternal parts, although they did not quantify the result. Contrary 
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Table 2. Continued. 









CN2 


C4668 


Fe4383 


Mgb 


Ca4227 


Galaxy 


grad(log age) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 




±cr 


±o- 


±0- 


±0- 


±(T 


±cr 


±0- 


NGC 3379 


0.1055 


-0.2279 


-0.515 


-0.419 


-0.239 


-1.069 


-0.262 




0.0347 


0.0528 


0.066 


0.051 


0.072 


0.128 


0.098 


NGC 3605 


0.0574 


-0.3890 


-0.050 


-0.385 


-0.869 


-0.452 


0.235 




0.1046 


0.1486 


0.218 


0.191 


0.278 


0.288 


0.348 


NGC 3608 


0.0595 


-0.3767 


-0.906 


-0.575 


-0.145 


-0.706 


0.104 




0.1043 


0.1180 


0.216 


0.209 


0.230 


0.306 


0.334 


NGC 3641 


0.0067 


-0.3821 


-1.608 


-0.547 


-0.209 


-0.856 


-0.261 




0.2529 


0.2634 


0.501 


0.359 


0.491 


0.636 


0.731 


NGC 3665 


0.0616 


-0.1746 


-0.092 


-0.383 


-0.105 


-0.086 


0.284 




0.0716 


0.0941 


0.166 


0.095 


0.163 


0.215 


0.285 


NGC 3818 


0.1865 


-0.3875 


-1.089 


-0.874 


-0.295 


-0.840 


-0.216 




0.0728 


0.1087 


0.152 


0.106 


0.147 


0.172 


0.203 


NGC 4261 


0.1011 


-0.2411 


-0.538 


-0.350 


-0.171 


-0.335 


-0.119 




0.0453 


0.0597 


0.075 


0.068 


0.091 


0.132 


0.184 


NGC 4278 


0.1895 


-0.4505 


-0.981 


-0.623 


-0.451 


-1.034 


-0.397 




0.0408 


0.0647 


0.081 


0.056 


0.095 


0.121 


0.132 


NGC 4365 


0.1087 


-0.2385 


-0.567 


-0.459 


-0.305 


-0.586 


-0.247 




0.0378 


0.0542 


0.069 


0.050 


0.077 


0.091 


0.121 


NGC 4374 


-0.1112 


-0.0716 


-0.452 


-0.258 


-0.000 


-0.239 


0.154 




0.0421 


0.0668 


0.109 


0.077 


0.092 


0.126 


0.158 


NGC 4415 


0.1048 


-0.0759 


0.074 


-0.038 


-0.071 


-0.428 


-0.581 




0.0625 


0.0658 


0.136 


0.080 


0.104 


0.243 


0.183 


NGC 4431 


0.1725 


-0.1285 


-0.610 


-0.216 


-0.260 


-0.409 


-0.548 




0.0928 


0.1048 


0.189 


0.145 


0.194 


0.219 


0.423 


NGC 4464 


-0.1041 


-0.0150 


-0.550 


-0.206 


-0.009 


-0.374 


0.213 




0.0796 


0.0931 


0.169 


0.141 


0.141 


0.178 


0.226 


NGC 4467 


0.0916 


-0.5234 


-0.798 


-1.202 


-0.719 


-0.669 


-0.303 




0.2661 


0.2767 


0.492 


0.504 


0.526 


0.782 


0.719 


NGC 4472 


0.1062 


-0.3919 


-0.543 


-0.472 


-0.220 


-0.331 


-0.506 




0.0252 


0.0460 


0.047 


0.035 


0.056 


0.072 


0.091 


NGC 4478 


0.1997 


-0.2858 


-0.508 


-0.428 


-0.391 


-0.393 


-0.454 




0.0508 


0.0667 


0.087 


0.056 


0.090 


0.122 


0.135 


NGC 4486B 


0.4092 


-0.3663 


-1.059 


-0.891 


-0.362 


-1.184 


-0.337 




0.2165 


0.2441 


0.260 


0.200 


0.483 


0.442 


0.576 


NGC 4489 


0.3077 


-0.4354 


-0.259 


-0.700 


-0.752 


-0.359 


-0.483 




0.1042 


0.1423 


0.228 


0.164 


0.260 


0.228 


0.492 


NGC 4552 


0.0526 


-0.2075 


-0.782 


-0.584 


-0.158 


-0.546 


-0.091 




0.0674 


0.0939 


0.127 


0.076 


0.151 


0.144 


0.287 



to these results, Mehlert et al. (2003) found, in a sample of 
galaxies from the Coma cluster, a mean age gradient compati- 
ble with being null. The discrepancy between the above results 
and those of Mehlert et al. (2003) may be due to the different 
environments from which the samples were drawn. We analyse 
the gradients for the HDEGs in Sec.|Sl 

Table 13 also shows the values of the rms dispersion about 
the mean value {cr) and the scatter expected from the errors 
(o"exp)- We performed a x^ test to compare these values and 
the results are shown in the 8* column of the table. There is 
a real scatter in the age gradients between galaxies that cannot 
be explained by the errors, consistent with the earlier results of 
G90. We have explored possible correlations between the age 
gradients and other parameters of the galaxies, including (but 
not limited to) the central velocity dispersion, but we have not 
found any significant correlations. 



4.2. Metallicity gradients 

The mean metallicity gradient in our sample of LDEGs is 
A[M/H]/Alogr=-0.206 + 0.019. Eariier studies have derived 
comparable gradients using, primarily (Fe) and Mg indices. For 
example. Couture & Hardy (1988), using the Mg2 index, found 
a variation in metallicity with radius of A[M/H]/ log r = -0.25. 
G90 measured mean gradients in their sample of early-type 
galaxies of A[M/H]/ log r = -0.23 ± 0.09 and A[M/H]/ log r = 
-0.22+0.10, using the Mg2 and (Fe) indices, respectively (later 
confirmed by Davies et al. 1993). To transform the line-strength 
gradients into metallicity gradients, each of these studies used 
the Mould (1978), Burstein (1979), and Faber et al. (1985) cal- 
ibrations. In all cases, they assumed a null age gradient. Fisher 
et al. (1995), through the comparison of the gradients of Mg2, 
H/3, and (Fe), with the stellar population models of Worthey 
(1994), obtained a metalhcity gradient of A[M/H]/logr = 



p. Sanchez-Blazquez et al.: Line-strength gradients in early-type galaxies 



11 



Table 2. Continued. 









CNj 


C4668 


Fe4383 


Mgb 


Ca4227 


Galaxy 


grad(age) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 




+o- 


±(r 


±(T 


±cr 


±cr 


+o- 


±<T 


NGC 4564 


0.1010 


-0.3264 


-0.842 


-0.740 


0.050 


-0.520 


-0.613 




0.0453 


0.0635 


0.107 


0.059 


0.080 


0.109 


0.131 


NGC 4594 


0.0938 


-0.3880 


-1.481 


-0.929 


-0.332 


-1.162 


-0.281 




0.1078 


0.1491 


0.272 


0.142 


0.335 


0.353 


0.496 


NGC 4621 


0.1645 


-0.4335 


-0.770 


-0.702 


-0.332 


-0.594 


-0.274 




0.0461 


0.0718 


0.097 


0.067 


0.114 


0.143 


0.184 


NGC 4636 


0.2275 


-0.3357 


-0.493 


-0.446 


-0.335 


-0.768 


-0.215 




0.0766 


0.0981 


0.178 


0.120 


0.158 


0.182 


0.289 


NGC 4649 


0.0305 


-0.2151 


-0.398 


-0.465 


-0.095 


-0.364 


-0.215 




0.0624 


0.0845 


0.159 


0.105 


0.142 


0.137 


0.261 


NGC 4673 


-0.1593 


-0.1049 


-0.231 


-0.540 


0.252 


-0.580 


-0.415 




0.0960 


0.1239 


0.176 


0.175 


0.213 


0.249 


0.339 


NGC 4692 


0.1630 


-0.3817 


-0.920 


-0.737 


-0.286 


-0.802 


-0.490 




0.1103 


0.1285 


0.233 


0.204 


0.251 


0.305 


0.349 


NGC 4697 


0.0425 


-0.1537 


-0.706 


-0.502 


-0.131 


-0.403 


-0.113 




0.0248 


0.0423 


0.046 


0.040 


0.055 


0.071 


0.086 


NGC 4742 


0.7591 


-0.6467 


-0.655 


-0.771 


-0.840 


-1.174 


-1.053 




0.0804 


0.0949 


0.129 


0.098 


0.143 


0.182 


0.207 


NGC 4839 


0.2800 


-0.3535 


-0.440 


-0.460 


-0.833 


-0.898 


0.056 




0.1708 


0.1540 


0.278 


0.220 


0.488 


0.647 


0.693 


NGC 4842A 


-0.2147 


0.0637 


-1.015 


0.098 


-0.324 


1.459 


1.735 




0.2617 


0.2364 


0.721 


0.295 


0.619 


0.676 


0.709 


NGC 4864 


0.4900 


-0.5792 


-0.479 


-0.697 


-1.175 


-0.998 


-0.692 




0.1762 


0.1618 


0.251 


0.244 


0.722 


0.812 


0.563 


NGC 4865 


0.0672 


-0.2587 


-0.410 


-0.359 


-0.143 


-0.532 


0.106 




0.1232 


0.1416 


0.370 


0.175 


0.261 


0.339 


0.368 


NGC 4874 


0.5845 


-0.7181 


-0.988 


-1.078 


-0.866 


-1.955 


-1.657 




0.1029 


0.1382 


0.225 


0.136 


0.265 


0.316 


0.411 


NGC 4875 


0.1777 


-0.3231 


-0.704 


-0.869 


-0.248 


-0.152 


-0.582 




0.1390 


0.2141 


0.852 


0.588 


1.212 


1.398 


1.876 


NGC 4889 


0.2578 


-0.4363 


-0.793 


-0.665 


-0.534 


0.052 


0.136 




0.0711 


0.1120 


0.182 


0.109 


0.163 


0.237 


0.307 


NGC 4908 


0.0131 


-0.2354 


-0.251 


-0.629 


-0.322 


-0.864 


0.613 




0.2231 


0.2812 


0.419 


0.604 


0.618 


0.549 


0.794 


NGC 5638 


0.0766 


-0.3273 


-0.774 


-0.515 


-0.140 


-0.614 


0.038 




0.0524 


0.0673 


0.097 


0.061 


0.103 


0.147 


0.187 


NGC 5796 


0.1234 


-0.2664 


-0.340 


-0.400 


-0.120 


-0.608 


-0.371 




0.0744 


0.0938 


0.133 


0.091 


0.148 


0.208 


0.216 



-0.25 ±0.1. Finally, Mehlert et al. (2003), using the population 
synthesis models of Thomas et al. (2003), obtained a metallic- 
ity gradient of A [M/H]/ log r = -0.10+0.12, significantly more 
shallow than previous studies had found. The values derived in 
the current study are consistent with the extant literature and 
suggest that early-type galaxies exhibit a reduction in metallic- 
ity of ~40% per decade of variation in radius. This value is con- 
siderably flatter than the values predicted by the dissipative col- 
lapse models. For example, Larson's hydrodynamical simula- 
tions gave AlogZ/Alog r 0.35 (Larson 1974a,b), and -1.0 

(Larson 1975) and Carlberg's (1984) N-body simulations gave 

A log Z/ A log r 0.5. It should be stressed though that these 

particular models leave significant room for improvement, due 
to the absence of essential physical processes, including star 
formation, thermal feedback from supernovae, and chemical 
enrichment. 



Finally, we compared the dispersion about the mean val- 
ues icr - 0.151) and the scatter expected from the errors 
(cexp=0.075) by performing a^^ test. The hypothesis "cTexp - 
cr" can be rejected at a low significance level (see 8* column 
of Table |3}. To investigate the causes of the variation in the 
metallicity gradients among galaxies, in Section|6l we explore 
putative correlations of the gradients with other physical galaxy 
parameters. 

5. Relative radial abundance ratios 

It is well known that massive early-type galaxies have abun- 
dance patterns that do not match that of the solar neighbour- 
hood (see Paper II, and references therein). Although the influ- 
ence of non-scaled solar abundance ratio patterns complicates 
the estimation of mean ages from integrated light, it does pro- 
vide an important clue as to the formation and chemical en- 
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Table 2. Continued. 









CNj 


C4668 


Fe4383 


Mgb 


Ca4227 


Galaxy 


grad(age) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 




±cr 


±(r 


±a- 


±cr 


±cr 


±o- 


±<T 


NGC5812 


0.1436 


-0.3077 


-0.795 


-0.567 


-0.355 


-0.480 


-0.434 




0.0645 


0.0881 


0.117 


0.079 


0.127 


0.172 


0.214 


NGC5813 


0.0150 


-0.0826 


-0.747 


-0.272 


-0.266 


-0.335 


-0.011 




0.0527 


0.0634 


0.105 


0.078 


0.113 


0.156 


0.202 


NGC5831 


0.0299 


-0.2412 


-0.827 


-0.608 


-0.417 


-0.553 


-0.401 




0.0682 


0.0980 


0.127 


0.123 


0.142 


0.149 


0.230 


NGC 5845 


0.0536 


-0.2156 


-0.591 


-0.482 


-0.540 


-0.337 


-0.304 




0.1005 


0.1117 


0.244 


0.150 


0.180 


0.216 


0.287 


NGC 5846 


0.1258 


-0.2327 


-0.432 


-0.363 


-0.786 


-0.389 


0.369 




0.0554 


0.0764 


0.130 


0.096 


0.165 


0.161 


0.236 


NGC 5846A 


0.1829 


-0.2979 


-0.570 


-0.848 


-0.917 


-1.069 


-0.451 




0.1110 


0.1620 


0.265 


0.156 


0.305 


0.278 


0.488 


NGC 6127 


0.2307 


-0.4596 


-0.879 


-0.667 


-0.261 


-0.500 


-0.258 




0.0855 


0.1293 


0.240 


0.118 


0.179 


0.375 


0.344 


NGC 6166 


-0.1318 


-0.0738 


-0.131 


-0.346 


-0.284 


-0.858 


-0.329 




0.0920 


0.1051 


0.137 


0.178 


0.228 


0.310 


0.271 


NGC 6411 


0.0348 


-0.1593 


-0.477 


-0.475 


-0.255 


-0.272 


-0.166 




0.0528 


0.0745 


0.107 


0.098 


0.122 


0.154 


0.181 


NGC 6482 


0.0031 


-0.1937 


-0.472 


-0.263 


-0.214 


-0.538 


-0.318 




0.0877 


0.0951 


0.170 


0.117 


0.185 


0.214 


0.321 


NGC 6702 


0.3652 


-0.5531 


-0.522 


-0.666 


-0.353 


-0.433 


-0.168 




0.0854 


0.1245 


0.204 


0.140 


0.150 


0.355 


0.343 


NGC 6703 


0.0898 


-0.2319 


-0.526 


-0.510 


-0.186 


-0.391 


0.112 




0.0423 


0.0684 


0.102 


0.059 


0.089 


0.093 


0.132 


NGC 7052 


-0.0011 


-0.1360 


-0.215 


-0.240 


-0.018 


-0.250 


0.309 




0.0596 


0.0763 


0.112 


0.105 


0.123 


0.153 


0.196 


NGC 7332 


-0.0207 


-0.0314 


-0.698 


-0.724 


-0.180 


-0.162 


0.217 




0.0460 


0.0768 


0.096 


0.075 


0.110 


0.113 


0.157 


IC767 


0.3838 


-1.1858 


-0.625 


-1.015 


-1.211 


-0.642 


-1.250 




0.1290 


0.2090 


0.459 


0.579 


1.359 


0.619 


0.683 


IC794 


0.2383 


-0.5238 


-0.230 


-0.831 


-0.580 


-0.314 


0.112 




0.1220 


0.1673 


0.229 


0.240 


0.236 


0.293 


0.472 


IC832 


-0.1432 


0.0902 


-0.329 


0.191 


0.242 


-0.073 


0.169 




0.1404 


0.1483 


0.291 


0.182 


0.330 


0.351 


0.710 



richment histories of galaxies. In this Section, we explore the 
presence of relative abundance ratios gradients in our sample 
of LDEGs To do so, we compare the metallicity gradients ob- 
tained with different indices (combined with HyS) and make the 
assumption that the differences in the derived values are due 
to the different sensitivity of the Lick indices to changes in the 
chemical composition. While avoiding a detailed quantitative 
analysis (due to the difficulty to perform this kind of analysis 
with the current stellar population models, see Paper II for de- 
tails), the qualitative trends presented below do offer invaluable 
information concerning the formation of these galaxies and the 
timescales for star formation therein. 

Fig.0shows three Aindex-Aindex diagrams in which gra- 
dients of CN2, Mgb and C4668 are compared against the gra- 
dients of Fe4383. The lines corresponding to Aage=0 and 
A[M/H]=0 are over-plotted. Due to the low sensitivity to age 
of all these indices, the diagrams appear highly degenerate - 
i.e., the iso-age and iso-metallicity lines are almost parallel. If 
the galaxies did not have a gradient in the relative abundances, 
we would expect to find all the points distributed along these 



lines. However, they appear to be systematically shifted toward 
the left of the diagram, although the magnitude of these shifts 
is very different in the three plots. In the first panel (AMgb- 
AFe4383), the metallicity gradients measured with Mgb are 
slightly steeper than the metallicity gradients inferred from 
Fe4383, although the differences are admittedly small. In the 
second panel (ACN2-AFe4383) the differences are more evi- 
dent. The metallicity gradients obtained with CN2 are clearly 
steeper (in absolute value) than the ones inferred from Fe4383. 
The last panel (AC4668-AFe4383) is an intermediate case be- 
tween the first two. The metallicity gradient obtained with 
C4668 is slightly steeper than the one estimated from Fe4383. 
Although the trend is not as evident as in the second panel, the 
points are visibly shifted toward the left of the constant age and 
metallicity lines. 

To analyse these differences in more detail, we calculated 
the metallicity gradients in several index-index diagrams using 
the indices CN2, C4668, Fe4383, Mgb and Ca4227, combined 
with Hy6. In Fig.|8] we compare these gradients with the values 
obtained from ten different indicators, as described in Sec. |3l 
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Table 2. Continued. 









CN2 


C4668 


Fe4383 


Mgb 


Ca4227 


Galaxy 


grad(age) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 


grad([M/H]) 




±(r 


±cr 


±0- 


±0- 


±0- 


±(r 


±0- 


IC 3957 


0.2298 


-0.4788 


-0.890 


-0.940 


-0.202 


-0.082 


-2.087 




0.2385 


0.2636 


0.769 


0.287 


0.507 


1.048 


0.612 


IC 3959 


0.0731 


-0.1642 


-1.038 


-0.077 


0.091 


-0.545 


-0.565 




0.1344 


0.1655 


0.214 


0.260 


0.349 


0.316 


0.510 


IC 3963 


0.3360 


-0.5703 


-0.424 


-0.864 


-0.550 


-0.765 


-1.322 




0.1340 


0.1788 


0.300 


0.217 


0.387 


0.436 


0.595 


IC 3973 


0.4122 


-0.3849 


-0.692 


-0.435 


-0.532 


-1.091 


-0.595 




0.1776 


0.1643 


0.374 


0.285 


0.483 


0.735 


0.820 


IC 4042 


-0.1029 


-0.3628 


-0.161 


-0.794 


0.919 


-0.275 


0.785 




0.4661 


0.6512 


0.374 


0.285 


0.483 


0.735 


0.820 


IC4051 


-0.0582 


-0.3693 


-0.873 


-0.759 


-0.224 


-0.365 


0.251 




0.0852 


0.1261 


0.165 


0.208 


0.184 


0.234 


0.411 


CGCG 159-41 


0.1629 


-0.0679 


-0.643 


0.156 


0.055 


-2.702 


0.558 




0.5742 


0.6968 


1.401 


1.034 


1.372 


1.443 


1.416 


CGCG 159-43 


0.5593 


-0.6726 


-0.902 


-0.751 


-0.480 


-1.305 


-0.820 




0.1614 


0.1936 


0.281 


0.355 


0.385 


0.386 


0.618 


CGCG 159-83 


-0.4099 


0.0939 


-0.902 


-0.751 


-0.480 


-1.305 


-0.820 




0.1310 


0.1566 


0.281 


0.355 


0.385 


0.386 


0.618 


CGCG 159-89 


-0.0815 


-0.3979 


-0.994 


-0.743 


0.570 


-0.048 


0.831 




0.1410 


0.1969 


0.229 


0.134 


0.351 


0.484 


0.513 



Table 3. Mean gradients of age, [M/H], and metallicity, derived with different indices for LDEGs. cr: typical deviation; A^: number 

of galaxies; A^eff: effective number of points, A^eff = Zi(l/c^) / 2 (1/c^)^ h t'- f-statistic to check the hypothesis "mean^O"; 
o"exp: typical deviation expected from errors; a: level of significance to reject the hypothesis "cr - o-gxp". The last column contains 
the f-statistic used to test the hypothesis "mean grad[M/H] = mean grad[M/H]" (with different indices). 





mean 


(T 


A' 


A'eff 


t 


Cexp 


a 


t 


grad(log age) 


0.082 ± 0.015 


0.117 


61 


19.62 


3.1 


0.032 


1.63E-13 




grad [M/H] 


-0.206 ± 0.019 


0.151 


61 


17.25 


5.7 


0.075 


4.3E-22 




grad [M/H] (CNj-H/?) 


-0.582 + 0.032 


0.249 


61 


23.14 


11.2 


0.111 


2.6E-31 


5.94 


grad [M/H] (C4668-H/3) 


-0.459 ± 0.028 


0.219 


61 


21.78 


9.8 


0.078 


0.000 


4.26 


grad [M/H] (Fe4383-H0) 


-0.197 ±0.025 


0.197 


61 


15.31 


3.9 


0.109 


1.2E-14 


0.14 


grad [M/H] (Mgb-H/3) 


-0.407 ± 0.039 


0.302 


61 


14.53 


5.1 


0.138 


5.3E-31 


2.30 


grad [M/H] (Ca4227-H/3) 


-0.238 ± 0.034 


0.269 


61 


19.18 


3.4 


0.178 


2.9E-7 


0.45 



Clearly, while the values obtained with Ca4227 and Fe4383 are 
compatible with the gradients obtained in Sec.|3l the gradients 
calculated with CN2, C4668 and, perhaps, Mgb are steeper We 
quantified these difference by measuring the mean gradient us- 
ing the various diagrams, the results of which are summarised 
in Table 13 The final column of the table shows the probability 
that the mean gradients calculated with the various indicators 
are the same as the gradients calculated Section |3 The gradi- 
ents calculated using CN2 and C4668 are almost twice as steep 
as the gradients calculated with the average often indicators. In 
the case of Mgb, despite the mean gradient being almost twice 
as large steep as the mean metallicity gradient, the result is less 
significant. In that case, the initial hypothesis can be rejected 
with a significance level lower than 0.025. 

We could argue that the differences in the metallicities de- 
rived with different indicators are due to variations in the chem- 
ical composition as a function of galactocentric distance, as the 
sensitivity of different Lick indices to variations on the abun- 
dance of different chemical elements is not the same (Tripicco 



& Bell 1995; Korn et al. 2004). But the Lick indices show a de- 
pendence with gravity that, while not large (Gorgas et al. 1993; 
Worthey et al. 1994), is also not null. As the behaviour of each 
line index under variations in the underlying IMF is not iden- 
tical, a change in the ratio of dwarf-to-giant stars with galac- 
tocentric radius could result in the metallicity gradient inferred 
from different indicators also being different. 

Cenarro et al. (2003) calibrated the CaT* (calcium triplet) 
index in the near-infrared and found a relationship between 
[Fe/H], velocity dispersion, and the slope of the IMF. One of 
the projections of this relationship is 



^i = 2.41 + 2.78[M/H] - 3.79[M/H]^ 



(12) 



where n represents the slope of the IMF. Using this equation we 
determined the slope of the IMF from the average metallicities 
in the centre, and at a distance of one effective radius, finding 
the following values: 

- /^ = 2.56 in the central regions, 
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Fig. 7. Aindex-Aindex diagrams in which we compare the Fe4383 gradients with the gradients of Mgb, CN2 and C4668. Dashed 
lines indicate the expected gradients if the only parameter varying with radius is the age (assuming an invariant solar metallicity), 
while solid lines show the expected trends if the only parameter changing with radius is the metallicity (assuming a constant age 
of 10 Gyr). 
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Fig. 8. Comparison of the metallicity gradient measured with different indicators combined with HyS, and the metallicity gradient 
calculated combining ten different indicators (see text for details). 



- H - 1.67 at a distance of one effective radius from the 
galaxy centre. 

To investigate the variation in the line-strength indices that 
such a variation in the IMF, plus an average variation in metal- 
licity as obtained in Section |3 would produce, we parame- 
terised the indices (using the V06 models) as a function of 
metalUcity and //, obtaining: 

CN#.0374 - 0.0063a/ + 0.1350[M/H] + 0.0495[M/H]2 

Ca422a.3745 + 0.1600yU + 0.9547[M/H] + 0.1566[M/H]2 

Fe438S.5963 + 0.0467;U + 4.2390[M/H] -1- 0.8785[M/H]2 

C466S.5930 - 0.2943// -1- 6.4247[M/H] + 1.8047[M/H]2 



MgS.7700 + 0.1303ju + 2.2567[M/H] + 2.2567[M/H]2 



A variation in the slope of the IMF, together with a varia- 
tion in metallicity equal to the mean metallicity gradient cal- 
culated in Section 12 would produce the following differences 
in the selected indices: ACN2 = 0.0278 mag, ACa4227= 
0.0390 mag, AFe4383= 0.0264 mag, AC4668=0.0182 mag, 
and AMgb=0.0270 mag. The mean observed variations for the 
LDEGs are: ACNzobs = -0.0662 + 0.0096 mag, ACa4227obs = 
-0.0157 ± 0.0097 mag, AFe4383obs = -0.0174 + 0.0064 mag, 
AC4668ob,s = -0.0314 + 0.0044 mag, and AMgbobs = 
-0.0301 + 0.0062 mag. As can be seen, a variation of the IMF 
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of the form predicted by Cenarro et al. (2003) cannot be re- 
sponsible for the different metallicity gradients obtained with 
the various indicators. In fact, this variation of the IMF would 
produce positive gradients in the selected indices. In this pa- 
rameterisation we have not included the effect of age, but the 
age gradients are not very strong in our sample (see Section|3 
and the differences in the sensitivities of the analysed indices 
to this parameter are not enough to produce the observed dif- 
ferences. 

More likely, the differences in the metallicity gradients ob- 
tained with the different indices are due to changes in the rela- 
tive abundance gradients of distinct elements. However, it is 
difficult to quantify the strength of those gradients with the 
present knowledge of line formation and the inherent limita- 
tions of stellar population synthesis models 

6. Correlation of the gradients with other 
parameters 

There are several physical processes that can produce a metal- 
licity gradient in early-type galaxies. The relationship between 
these gradients and other global properties of the galaxies af- 
ford an opportunity to discriminate between the competing pro- 
cesses. Previous attempts have searched for correlations using 
colours and empirical line-strength indices; in what follows, we 
study the correlation between the gradients of the derived sim- 
ple stellar population (SSP) parameters and the global prop- 
erties of the host galaxies. In this section, we only use the 
subsample of elliptical galaxies, excluding the lenticulars (SO) 
from the analysis. The reason for doing this is that the cor- 
relations between the gradients and other parameters are pre- 
dicted by an specific mechanism of galaxy formation, and, al- 
though SO and E galaxies seem to follow the same relations 
between the central properties and other parameters, very dif- 
ferent mechanism have been proposed for their formation. As 
the sample of SO is small, we could not perform a comparative 
study using only SO galaxies. 

6.1. Correlation of the metallicity gradient with the 
velocity dispersion gradient 

As noted in Section 1, Franx & Illingworth (1990) found a cor- 
relation between colour gradient and the gradient of escape 
velocity, which they interpreted as a correlation between the 
local metallicity and the local potential well depth in a sam- 
ple of early-type galaxies. They suggested that such a corre- 
lation was consistent with a galactic wind origin, and incon- 
sistent with a dissipative inward flow origin. Others, including 
Davies et al. (1993), Carollo & Danziger (1994), and Mehlert et 
al. (2003) find similar relations between the gradients of some 
line-strength indices and the gradient of the potential well.^. 
We now examine the likelihood for the existence of this corre- 



■^ Davies et al. (2003) calculated the local escape velocity using 
the surface brightness and the kinematics of the galaxies (J0rgensen, 
Franx & Kaergaard 1992). Carollo & Danziger (1994) obtained this 
parameter with asymmetric dynamical models which depended on the 
total energy and the angular momentum along the symmetry axis. 



lation within our sample of LDEGs, using the infeiTed metal- 
licity gradients and empirical velocity dispersion gradients - 
Fig-Elillustrates the observed trends. 

To study the degree of coiTelation we have performed two 
different tests: a non-parametric Spearman rank order test and 
a f-test. The Spearman test does not take into account the er- 
rors in the individual points, while the f-test has the limitation 
of considering just a linear relation and it assumes Gaussian 
probability distributions. A good estimate of the degree of cor- 
relation between both magnitudes can be obtained by studying 
the results of both tests. In the f-test, we check the hypothe- 
sis "^7 - 0", where b is the slope of the linear fit to the data. 
The errors in the slope were calculated through Monte Carlo 
simulations, where each point was perturbed in both axes as- 
suming a Gaussian distribution with standard deviation equal 
to the errors. For a significance level of a = 0.05, a t value 
higher than 1.96 indicates that a correlation exists. The results 
of both tests are shown in Table|3 In the second row of the ta- 
ble (and left top panel in Fig.|9j we show the correlation when 
the metallicity is measured with the ten indicators described in 
Sec. 13 We do not find any congelation between the metallic- 
ity gradients and the gradient of velocity dispersion. The other 
panels in the figure show the metallicity gradients obtained in 
different Aindex-Aindex diagrams in which we combined H/3 
with various metallicity indicators (as indicated between the 
brackets). Surprisingly, we find a marginally significant corre- 
lation between the metallicity gradients and grad cr when the 
metallicity gradient is measured using Mgb and CN2 indices, 
while we do not find any correlation when the metallicity gra- 
dient is obtained using the other indices. If the correlation be- 
tween these parameters is indicative of the importance of galac- 
tic winds during the evollution of early-type galaxies, the fact 
that the correlation only exists when the metallicity gradient 
is calculated with some specific indices may indicate that this 
process affects some chemical species more than others. This 
could happen if, for example, the relative abundance patterns 
were different when the galactic winds occurred. 

The dispersion in the relations is, in any case, very large. 
This large scatter may be a consequence of using the velocity 
dispersion as an indicator of the local potential and/or indica- 
tive that other processes are driving the variation of the local 
metal content with radius. Davies et al. (1993), in fact, suggest 
that the local velocity dispersion is a poor indicator of the local 
escape velocity, due to the (complicating) presence of rotation 
and anisotropics. On the other hand, the existence of kinemat- 
ically decoupled cores in a large percentage of galaxies, the 
presence of shells, dust lanes, and the observation of interact- 
ing galaxies, seem to indicate that mergers and interactions are 
common processes in the lives of galaxies. If these interactions 
have some associated gas dissipation and/or star formation, it 
may produce dispersion in the correlation between local metal- 
licity and local escape velocity. In fact, other authors (Davies et 
al. 1993; Carollo & Danziger 1994) have found correlations be- 
tween the Mg2 and the velocity dispersion gradients using the 



Mehlert et al. (2003) used the gradient of velocity dispersion as a mea- 
sure of the local potential well depth. 
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Fig. 9. Relation between the metallicity gradients obtained in different index-index diagrams (as indicated in the brackets) and 
the velocity dispersion gradient for the sample of LDEGs. In the top left panel, the metallicity gradient has been obtained using 
a combination often indices, as described in Sec.|3l 

Table 4. Linear fits and probability of no correlation between the metallicity gradient and the gradient of velocity dispersion. A^: 
number of galaxies in the fit; b + a-{b): slope and error of the linear fit (the errors have been calculated by Monte Carlo simulations 
as described in the text); t: f-statistic to verify the hypothesis "Z? = 0" (a t value higher than 1.96 allows the hypothesis to be 
rejected, with a significance level lower than 0.05); Pnc: probability of no coiTelation in a non-parametric Spearman test. 





A' 


b + o-{b) 


t 


Pnc 


grad age 


54 


1.727 + 3.099 


0.56 


0.885 


grad [M/H] 


54 


0.175 ±0.269 


0.65 


0.238 


grad [M/H] (CN^-H/J) 


54 


1.134 ±0.555 


2.04 


0.002 


grad [M/H] (Ca4227-H/3) 


54 


0.449 ± 0.686 


0.65 


0.063 


grad [M/H] (Fe4383-H/3) 


54 


-0.392 ± 0.442 


0.89 


0.900 


grad [M/H] (C4668-Hy3) 


54 


0.451 ±0.356 


1.27 


0.062 


grad [M/H] (Mgb-HyS) 


54 


0.876 ± 0.463 


2.11 


0.001 



local escape velocity instead of the velocity dispersion, finding 
also a large scatter in the relations. 

We conclude that the local potential of early-type galaxies 
may play a role in defining the metallicity gradient. However, 
the large scatter in the derived coiTelation suggests that other 
processes may play a role in modulating the final metal content. 
These processes can be the consequence of differences in the 
merger histories of galaxies as proposed by Kobayashi (2004). 
The fact that the inferred metallicity gradients differ depending 
upon the indicator adopted also suggests that the gradients of 
the various chemical species have probably been formed (and 
modified) through different physical mechanisms. 

6.2. Correlation of the metallicity gradient with the 
central velocity dispersion 

Dissipative collapse models of galaxy formation predict strong 
coiTelations between the metallicity gradients and certain 
global parameters, such as the luminosity, mass, and central 
velocity dispersion (Larson 1974a; Carlberg 1984; Arimoto & 
Yoshii 1987; Kawata 1999; Chiosi & Carraro 2002; Kawata 



& Gibson 2003), in the sense that more massive galaxies 
should possess steeper metallicity gradients. Such a predic- 
tion is driven primarily by the adopted mass-dependent feed- 
back efficiency within the models. Conversely, galaxy forma- 
tion through hierarchical clustering of small sub-units does not 
necessarily lead to a clear prediction for any putative correla- 
tions between metallicity gradients and global galactic prop- 
erties - merger history can readily play a part in eroding any 
extant correlation (e.g. Kobayashi 2004). 

With the aim of exploring the possible correlation between 
metallicity gradient and galactic mass, we compare in Fig. 1101 
the metallicity gradients (obtained as described in Sec. O with 
the central velocity dispersion for the sample of LDEGs. The 
other panels in the figure show the relation between the metal- 
licity gradient and the central velocity dispersion when the 
metallicity gradient is measured using CN2, C4668, Fe4383, 
Mgb and Ca4227 combined with H/3. 

As in Section 6.1, we studied the degree of coiTelation 
through a f-test and a non-parametric Spearman test. The re- 
sults of these tests are shown in Table |5] We did not find any 
correlation between the metallicity gradients and the central ve- 
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Table 5. Correlation between the gradients of age and metal- 
licity and the central velocity dispersion. The brackets show 
the different indicators used to calculate the metallicity gradi- 
ents. The first and the second rows refer to the values calculated 
using the ten indicators described in Section |3l t: t-parameter 
calculated as the quotient of the slope and its error, computed 
through Monte Carlo simulations; Pnc: probability of no corre- 
lation according to a non-parametric Spearman test. 
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locity dispersion, regardless of the line-strength indices used to 
derive the former. 

In regards to putative correlations between line-strength in- 
dices and the mass or the velocity dispersion, Peletier (1989), 
Peletier et al. (1990), and G90, each found positive correlations 
between the Mg2 gradient and the central velocity dispersion, 
while Gonzalez (1993) and Davies & Sadler (1987) found the 
opposite trend. Conversely, Davidge (1991, 1992) and Davies 
et al. (1993) did not find any correlation between the gradi- 
ents of the indices Mg2, <Fe> or Hfi and the central velocity 
dispersion. In this work, we have doubled the sample size com- 
pared to any of the extant stuudies, and instead of analysing the 
colours or raw line-strength gradients, we have employed the 
inferred metallicity gradients. For the sample of LDEGs, we 
did not find any correlation between these parameters and the 
central velocity dispersion. 



6.3. Correlation of the gradients witli the age of the 
central regions 

In Section I6J1 we showed that, for the sample of LDEGs, the 
metallicity gradients tend to correlate with the velocity disper- 
sion gradient when they are measured with some indicators 
(in particular CN2 and Mgb) but that this correlation disap- 
pears when the metallicity gradient is measured with Fe4383, 
Ca4227, and C4668. 

In Section|3]we showed that the mean age gradient in our 
sample of LDEGs is not null, which might be indicative of 
the occurrence of star formation events in the centre of the 
galaxies. Furthermore, we showed in Paper II than in these 
events, the relative enrichment of some chemical species (the 
ones mainly released by low- and intermediate-mass stars, like 
Fe) might be more important than others (the ones mainly pro- 
duced in Type II supernovae, like Mg). With the aim of explor- 
ing if the occurrence of star formation processes in the centre 
of the galaxies may be responsible for the abundance gradi- 
ents in some chemical species, we now study the correlation 
between the metallicity gradients and the central ages obtained 
in Paper II. Fig.^Jshows these correlations for different mea- 
surements of the metallicity gradient. When the correlation is 
statistically significant, a linear fit taking into account the er- 
rors in the x- and y-directions is also plotted. In these cases, the 
probability of no correlation (Pnc) obtained in a non-parametric 
Spearman test, and the t parameter to verify the hypothesis 
b - (where b is the slope of the linear fit), are indicated 
in the panels. As can be seen in the different panels, while we 
do not find any significant correlation when the metallicity gra- 
dients are measured with CN2, Mgb, or Ca4227, we do find 
a correlation between the metallicity gradients and the central 
age of the galaxies when the former is measured using Fe4383 
and C4668. We also find a significant correlation between the 
metallicity gradient obtained with the ten indicators described 
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Fig. 11. Relation of the metallicity gradient with the central age of the galaxies. When the correlation is statistically significant 
a linear fit weighted by the errors in both axes is also shown. In these cases, the t parameter and the probability of no correla- 
tion (Pnc), as derived from a non-parametric Spearman rank test, are also indicated in the top-left corner of the panels. For a 
significance level of 0.05, a value of t larger than 1 .96 indicates the presence of correlation between both variables. 



in Section |3l and the central age of the galaxies. We interpret 
these correlations in Section 6.4. 



6.4. Interpretation of the Correlations for the LDEGs 

In the Sections 6.1-6.3, we studied the putative correlations 
between metallicity gradient and other global properties of 
LDEGs. Broadly speaking, we found two distinct characteris- 
tics of the gradients depending upon the specific indicator used 
to derive them: 



- Metallicity gradients measured with CN2 and Mgb. These 
correlate with the velocity dispersion gradient but they do 
not correlate with the central velocity dispersion or with 
the central age of the galaxies. This behaviour may indi- 
cate that the gradients of some elements (in particular Mg 
and N) are the result of processes such as galactic winds 
which can develop after an intense burst of star formation 
in the earliest phases of galaxy formation. The possible sub- 
sequent star formation in the galaxies has not modified sub- 
stantially the variation of these elements with radius. This 
is in agreement with the results of Papers I and II where we 
argued that the relative enrichment of Fe in the star forma- 
tion processes which occurred in the centre of the galaxies 
had to be more important than that of Mg. 

- Metallicity gradients measured with Fe4383 and C4668. 
These gradients do not correlate with the velocity disper- 
sion gradient but they do correlate with the central age of 
the galaxy. The correlation with the central age may in- 
dicate that episodes of star formation in the centre of the 
galaxies affects considerably the metallicity gradient for 
these chemical species, increasing the central metallicity. 



This view agrees with the results of Papers I and II. Note 
that in galaxies that have suffered star formation episodes 
in their centres, the [Mg/Fe] gradient flattens as a conse- 
quence of the Fe enrichment in the central parts, which 
is supported by the relation between [Mg/Fe] and age re- 
ported in Paper I. This is also in agreement with the result 
presented in Paper II in which we found the existence of an 
age-metallicity relation when the metallicity was measured 
with Fe4383, but not when this parameter was measured 
with Mgb. 

Conversely, in Sanchez-Blazquez et al. (2003) we reported 
differences in the N abundance between galaxies in differ- 
ent environments. If the star formation processes in the cen- 
tres of the galaxies do not affect substantially the gradients 
in this element, then the differences must be visible not only 
in the centres, but in the outer parts. 

We also note the behaviour of the metallicity gradient when 
inferred from Ca4227. Specifically, the gradient does not 
correlate with any of the obvious physical parameters - i.e., 
neither with the central age nor with the velocity dispersion 
gradient. 

7. Global stellar population parameters 

Many of the results described in Papers I and II could be ex- 
plained by assuming the presence of a small percentage of 
young stars in the centre of most of the galaxies, at least in 
the subsample of LDEGs. We have also argued that this could 
be the cause of the existence of non-null age gradients in our 
sample of galaxies (see Sec. |4}. If this point of view is cor- 
rect, we would expect that the relations defined between the 
global parameters of the galaxies (i.e. for the whole bodies of 
the galaxies) are different than the one derived for the central 
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regions. To investigate this possibility, we now compare for the 
subsample of LDEGs the relation of these global values and 
the velocity dispersion with the relations derived for the cen- 
tral values in Paper II. The global values can be obtained from 
the gradients, assuming a linear behaviour of the indices with 
radius and evaluating the integral 
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where I{r) - a + blog(r/re) represents the index at a distance 
r from the galaxy as derived from the gradient, and r^ is the 
effective radius. To solve this integral, we have assumed that 
the spatial profile of the galaxies can be approximated with a 
de Vaucouleurs law. 

Fig. ^] shows the relation between the global age and 
global metallicities (as derived from different indicators) and 
the central velocity dispersion. The solid line indicates a lin- 
ear fit to the data weighted with the errors in both parameters. 
In order to compare with the relations for the central regions, 
we have also plotted the linear fits obtained for the central val- 
ues derived in Paper II (dashed lines). Table |6l summarises the 
parameters of the fits. We carried out a f-test to verify the exis- 
tence of significant differences between the slopes of the trends 
defined by the central and the global SSP-parameters. A value 
of t higher than 1.96 indicates that there exist differences with 
a significance level lower than 0.05. 

The only case in which the slopes defined by the central 
and the global values are significantly different is in the rela- 
tion between the age and the central velocity dispersion. While 
the central age shows a significant correlation with the central 
velocity dispersion (see Paper II), the global age as derived in 
this section does not correlate with this parameter. This result 
favours the idea that the young ages that we found in Paper II 
in a large percentage of galaxies are due to a minor percentage 
of stars formed in the centre of the galaxies at a later epoch 
than the bulk of the stars, as suggested by Trager et al. (2000b). 
It also supports the suggestion that this minor percentage of 
young stars in the centre of the galaxies is responsible for the 
age gradients reported in Section|3 

On the other hand, the slopes of the relations between the 
metallicity and the central velocity dispersion do not show a 
statistically significant variation between the central and the 
global values, although there is a tendency for the relations to 
be steeper for the global values. There are also differences in 
the zero point as a consequence of the existence of gradients, 
in disagreement with the earlier claims of Gonzalez & Gorgas 
(1996). These authors found that the relation Mg2 - cr is flat- 
ter at one effective radii than in the center, concluding that the 
mass-metallicity relation was much flatter at one effective radii 
than in the central parts of the galaxies. In light of our present 
series of papers, the results of Gonzalez & Gorgas could be 
explained if the differences between the central and global re- 
lation of the age with cr was causing these differences. 

The lack of variation in the slope of the relation between 
the central and the global metallicities measured with CN2 and 
Mgb was expected, since we find a flat relation between the 
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Fig. 12. Relation between the global age and metallicity and the 
central velocity dispersion for galaxies in low-density environ- 
ments. The solid line represents a linear fit weighted with the 
errors for the global data. The dashed line indicates the fit ob- 
tained for the central values reported in Paper II. Open circles 
represent E galaxies while the SO are represented by squares. 

strength of the gradients and the central velocity dispersion. 
In the case of the metallicity measured with Fe4383 we might 
have expected to see a steepening of the slope in the global re- 
lation compared with the central one. If the relative importance 
of the star formation processes has been higher in the smaller 
galaxies (as suggested by the age-cr relation), the metallicity 
inferred from Fe4383 in these galaxies should be also higher. 
The variation in the slope obtained when comparing the central 
and the global relations goes in this sense, but the differences 
are not statistically significant. We must note that the errors in 
the global measurements are higher than for the central values, 
which might explain the lack of statistical significance. 

8. Differences in the gradients as a function of 
environment 

If the environment in which galaxies reside has any influence 
over the timescales of the star formation, the number of inter- 
actions, or the dissipation of gas, we might expect to see an 
environmental dependence upon the inferred age and metallic- 
ity gradients. Fisher et al. (1995) analysed a sample of bright 
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Table 6. Parameters of the linear fits, weighted with errors, of the age and metallicities (measured in different index-index 
diagrams indicated between brackets) and the central velocity dispersion. The table shows the results for the central values (oc 
and be) and for the global values (Og and b^). The coefficients a and b represent the zero point and the slope of the linear fit, 
respectively. The last column shows the t parameter obtained in a f-test to check the hypothesis be = b^. A value of t higher than 
2.326 allows us to reject the hypothesis with a significance level lower than 0.01. 
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cluster ellipticals, concluding that their gradients and those of 
field ellipticals, were not significantly different. Tamura & Otha 
(2003) found, studying the photometry in the B and R bands as 
a function of radius for galaxies in Abell 2199, correlations 
between the colour gradients and some global properties of the 
galaxies, such as the luminosity and effective radius, which had 
not been found in studies of field galaxies. These authors, how- 
ever, only found these correlations amongst the most luminous 
galaxies (R < 15 mag and with an effective radius >3")- We 
now we examine the mean gradients of the SSP-parameters for 
the HDEGs and compare these values with the ones obtained 
for the LDEGs. 

8.1. Mean stellar population gradients in HDEGs 

The second column of Table lists the mean values of the 
SSP-parameter gradients obtained for the HDEGs in differ- 
ent index-index diagrams (indicated between brackets). As in 
Section 0] we quantified the probability that these values are 
different from zero by chance. The results of a f-test indi- 
cate that, while the metallicity gradients obtained in the CN2- 
Hy6 and C4668-HyS diagrams are not compatible with being 
null, the metallicity gradients obtained in the Fe4383-Hy6 and 
Ca4227-Hy6 diagrams are compatible with zero, within the er- 
rors. Furthermore, the mean age gradient for this subsample of 
galaxies is also compatible with zero. We now analyse these 
results separately: 

- Age gradient. If we assume that the most likely scenario 
to explain the age gradient in early-type galaxies is the oc- 
currence of star formation processes in the centres of these 
systems, the lack of a mean age gradient in the HDEGs in- 
dicates that these galaxies have undergone few episodes of 
star formation in recent times, when compared with field 
galaxies. If the star formation processes are triggered by 
the interactions between galaxies, the differences could be 
explained due to the lower probability of an interaction in 
the centre of the Coma cluster. This is in agreement with 
the conclusions of Papers I and II, where we showed that 
the central stellar populations of the Coma cluster galaxies 
possessed, on average, older ages than the LDEGs. 

- Metallicity gradients. The mean metallicity gradient ob- 
tained for this subsample of galaxies is A[M/H]/logr = 
-0.328 + 0.064, slightly steeper than the gradient obtained 
for the LDEGs. From her simulations, Kobayashi (2004) 



predicts a metallicity gradient for galaxies that have not 
suffered major mergers of AlogZ/Alogr ~ -0.3. These 
predictions are in agreement with the mean values obtained 
for the galaxies in the Coma cluster On the other hand, as 
for the LDEGs, the metallicity gradients are steeper when 
obtained with some indicators (CN2, C4668 and Mgb), al- 
though the statistical significance of the differences is much 
lower than in the case of the LDEGs (see final column of 
Table 0. This may also be due to the larger errors in the 
determination of the gradients for the galaxies in the Coma 
cluster. 

We have next checked to see if the scatter amongst the mean 
values was compatible with the dispersion expected by the er- 
rors. We carried out ax^ test of the hypothesis cr = cTexp (where 
cr is the observed scatter and o-f-^p the scatter expected from 
the errors), to see whether it could be rejected with a low sig- 
nificance level {a). Column 8 of Table shows the a values 
resulting from this test. We find a real scatter in both the age 
and the metallicity gradients which cannot be explained by the 
errors, in contrast with the findings of Mehlert et al. (2003). 
The difference may be due to the fact that the sample used by 
Mehlert et al. spans a more limited range in velocity dispersion 
(2.2 < logo- < 2.5). 

8.2. Correlation of the metallicity gradients in HDEGs 
with the central velocity dispersion 

After Kobayashi (2004), the absence of correlation between 
the metallicity gradient and the central velocity dispersion may 
be the consequence of differences in the merger history of the 
galaxies. If the star formation processes in recent epochs have 
been less frequent in the Coma cluster galaxies, one might ex- 
pect that any in situ correlation would have been impacted upon 
less substantially. To explore if this is the case, in FigureFT^we 
present the metallicity gradients measured in different Aindex- 
Aindex diagrams against the central velocity dispersion for the 
elliptical galaxies of the Coma cluster. The table at the bottom 
of the figure shows the results of a f-test and a non-parametric 
Spearman rank test to check the degree of correlation between 
both variables. In this case we find a correlation between the 
metallicity gradient and the central velocity dispersion, but 
only when the metalUcity gradient is measured with the Fe4383 
index. While the correlation is not statistically significant when 
derived from other indices, there does seem be a marginal trend 
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Table 7. Mean stellar population gradients and their errors for HDEGs. o". standard deviation about the mean; A^: number of 
galaxies averaged; A^eff : effective number of points; f: f-statistic to verify the hypothesis "mean=0"; cr^,^^: standard deviation 
expected from the errors; a: significance level to reject the hypothesis "cr^CTexp" in a. x^ test. Column 9 shows the f-statistic to 
confirm if the gradients obtained with individual indicators are the same as the metallicity gradients obtained with a combination 
of ten different indicators. As a reference, the 10* column shows the mean gradients obtained for the LDEGs, from Section]?] 
(see Tabled- The final column of the table shows the f-statistic, to verify if the mean gradients are the same for both samples of 
galaxies, LDEGs and HDEGs. A high value of f indicates significant differences. 
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in the sense that more massive galaxies also show a steeper gra- 
dient. To confirm this putative mass trend and further constrain 
galaxy formation models, higher signal-to-noise spectra must 
be obtained. 

9. Conclusions 

We have carried out a study of the gradients in 23 different 
spectral features for a sample of 82 early-type galaxies situ- 
ated in different environments. Our results can be summarised 
thusly: 

- Using the new synthesis stellar population models of V06 
we have derived age and metallicity gradients for all the 
galaxies in the sample. We have used a new method with 
employs ten different indicators in order to reduce the scat- 
ter due to random errors. 

- The mean age and metallicity gradients for the LDEGs are 
A[M/H]/logr = -0.205 + 0.075 and A log(age)/log r = 
0.082 + 0.032, respectively. The mean age gradient is 
steeper and the mean metallicity gradient flatter, than the 
predictions of dissipative collapse models. On the other 
hand, the dispersion amongst the mean values is larger than 
the dispersion expected by errors alone. 

- We have studied the metallicity gradients derived from dif- 
ferent indicators, obtaining steeper gradients when using 
CN2 and C4668 than with Fe4383 and Ca4227. Although it 
is beyond the scope of this paper to derive chemical abun- 
dances ratio gradients, we speculate that these aforemen- 
tioned differences are reflecting the existence of radial vari- 
ations in the relative abundances of some elements with re- 
specto to iron. 

- We have studied the relation between the age and metallic- 
ity gradients and the gradient of velocity dispersion, find- 
ing a correlation when the metallicity gradient is calculated 
with CN2 and Mgb. We do not, however, find a correlation 
when the metallicity gradient is derived from the Fe4383 
and Ca4227 indices. The fact that the gradients of CN2 
and Mgb correlate with the gradient of the velocity disper- 
sion may indicate that the gradients in some elements (e.g. 
N and Mg) were shaped early in the formation of galax- 



ies, when the galactic winds presumably dominate. On the 
other hand, the lack of correlation when using Fe4383 and 
Ca4227 to measure the metallicity can indicate that other 
process, such as secondary bursts of star formation, have 
had a stronger influence in the final shape of these gradi- 
ents. 

- We have not found any correlation between the metallicity 
gradients and the central velocity dispersion for the sample 
of galaxies in low-density environments. A strong correla- 
tion between the gradients and the mass of the galaxy is 
expected within dissipative collapse formation scenarios. 

- We have found a significant correlation between the metal- 
licity gradients and the central age for the LDEGs when the 
metallicity gradients are measured with Fe4383 and C4668, 
in the sense that galaxies with a younger central age also 
show a steeper metallicity gradient. 

- The results quoted above suggest that the gradients of 
different chemical species may have formed by different 
mechanisms. 

- The mean age and metallicity gradients for the galaxies in 
the Coma cluster are Alog(age)/logr = 0.027+0.138 and 
A[M/H]/logr = -0.328 + 0.157, respectively. The mean 
value of the age gradient is compatible, to within the er- 
rors, with zero. Both values are also compatible with the 
predictions of Kobayashi (2004) for those galaxies which 
have not undergone major mergers. However, the disper- 
sion amongst the mean values is higher than the scatter ex- 
pected by the errors, indicative of real differences between 
galaxies. 

- For galaxies in the Coma cluster, we have studied the cor- 
relation between the metallicity gradients and the central 
velocity dispersion finding a statistically significant cor- 
relation between both parameters when the metallicity is 
measured with Fe4383. For the other indices, the correla- 
tion is not statistically significant, but we observe a trend in 
the sense that more massive galaxies tend to have a steeper 
metallicity gradient. This trend is predicted by models of 
dissipative collapse, and was not found for the LDEGs. We 
should stress though that the quality of the data for the 
HDEGs is lower than that for the LDEGS; higher quality 
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Fig. 13. Metallicity gradients versus central velocity dispersion for HDEGs. The different panels show the metallicity gradients 
obtained in different index-index diagrams (indicated between the brackets). The first panel (top left) compares the metallicity 
gradients obtained as described in Sec.|3] In the case in which the correlation is statistically significant, an error weighted linear 
fit is also shown. The table at the bottom indicates the results of a f-test and a non-parametric Spearman rank test (see text for 
details). 



data for the Coma cluster galaxies are urgently needed. In 
any case, in this paper we have found systematic differences 
between the stellar populations of galaxies in different envi- 
ronments which both confirms and extends our conclusions 
from Papers I and II. 
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